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Foreword 



id , 



This Technical Specification has been produced by the 3 Generation Partnership Project (3GPP). 

The present document specifies the speech codec to be used for the GSM half rate channel for the digital cellular 
telecommunications system. The present document is part of a series covering the half rate speech traffic channels as 
described below: 

GSM 06.02 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Half rate speech 

processing functions". 

GSM 06.06 "Digital cellular telecommunications system (Phase 2+); Half rate speech; ANSTC code for the 

GSM half rate speech codec". 

GSM 06.07 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Test sequences for the 

GSM half rate speech codec". 

GSM 06.20 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Half rate speech 
transcoding". 

GSM 06.21 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Substitution and muting 

of lost frames for half rate speech traffic channels". 

GSM 06.22 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Comfort noise aspects 

for half rate speech traffic channels". 

GSM 06.41 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Discontinuous 

Transmission (DTX) for half rate speech traffic channels". 

GSM 06.42 "Digital cellular telecommunications system (Phase 2+); Half rate speech; Voice Activity Detector 
(VAD) for half rate speech traffic channels". 

The contents of the present document are subject to continuing work within the TSG and may change following formal 
TSG approval. Should the TSG modify the contents of the present document, it will be re-released by the TSG with an 
identifying change of release date and an increase in version number as follows: 

Version x.y.z 

where: 

x the first digit: 

1 presented to TSG for information; 

2 presented to TSG for approval; 

3 or greater indicates TSG approved document under change control. 

y the second digit is incremented for all changes of substance, i.e. technical enhancements, corrections, 
updates, etc. 

z the third digit is incremented when editorial only changes have been incorporated in the document. 



ETSI 



3GPP TS 46.020 version 7.0.0 Release 7 6 ETSI TS 146 020 V7.0.0 (2007-06) 



Scope 



The present document specifies the speech codec to be used for the GSM half rate channel. It also specifies the test 
methods to be used to verify that the codec implementation complies with the present document. 

The requirements are mandatory for the codec to be used either in GSM Mobile Stations (MS)s or Base Station Systems 
(BSS)s that utilize the half rate GSM speech traffic channel. 



References 



The following documents contain provisions which, through reference in this text, constitute provisions of the present 
document. 

• References are either specific (identified by date of publication, edition number, version number, etc.) or 
non-specific. 

• For a specific reference, subsequent revisions do not apply. 

• For a non-specific reference, the latest version applies. In the case of a reference to a 3GPP document (including a 
GSM document), a non-specific reference implicitly refers to the latest version of that document in the same 
Release as the present document. 

[1] GSM 06.02: "Digital cellular telecommunications system (Phase 2+); Half rate speech; Half rate 

speech processing functions". 

[2] GSM 06.06: "Digital cellular telecommunications system (Phase 2+); Half rate speech; ANSI-C 

code for the GSM half rate speech codec". 

[3] GSM 06.07: "Digital cellular telecommunications system (Phase 2+); Half rate speech; Test 

sequences for the GSM half rate speech codec". 



3 Definitions, symbols and abbreviations 

3.1 Definitions 

For the purposes of the present document, the following definitions apply: 

adaptive codebook: adaptive codebook is derived from the long term filter state. The lag value can be viewed as an 
index into the adaptive codebook. 

adaptive pitch prefilter: in the GSM half rate speech decoder, this filter is applied to the excitation signal to enhance 
the periodicity of the reconstructed speech. Note that this is done prior to the application of the short term filter. 

adaptive spectral postfilter: in the GSM half rate speech decoder, this filter is applied to the output of the short term 
filter to enhance the perceptual quality of the reconstructed speech. 

allowable lags: set of lag values which may be coded by the GSM half rate speech encoder and transmitted to the GSM 
half rate speech decoder. This set contains both integer and fractional values (see table 3). 

analysis window: for each frame, the short term filter coefficients are computed using the high pass filtered speech 
samples within the analysis window. The analysis window is 170 samples in length, and is centered about the last 100 
samples in the frame. 

basis vectors: set of M, Ml, or M2 vectors of length Ns used to generate the VSELP codebook vectors. These vectors 
are not necessarily orthogonal. 

closed loop lag search: process of determining the near optimal lag value from the weighted input speech and the long 
term filter state. 
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closed loop lag trajectory: for a given frame, the sequence of near optimal lag values whose elements correspond to 
each of the four subframes as determined by the closed loop lag search. 

codebook: set of vectors used in a vector quantizer. 

Codeword (OR Code): M, Ml, or M2 bit symbol indicating the vector to be selected from a VSELP codebook. 

Delta (LAG) code: four bit code indicating the change in lag value for a subframe relative to the previous subframe's 
coded lag. For frames in which the long term predictor is enabled (MODE 1, 2, or 3), the lag for subframe 1 is 
independently coded using eight bits, and delta codes are used for subframes 2, 3, and 4. 

direct form coefficients: one of the formats for storing the short term filter parameters. All filters which are used to 
modify speech samples use direct form coefficients. 

fractional lags: set of lag values having sub-sample resolution. Note that not every fractional lag value considered in 
the GSM half rate speech encoder is an allowable lag value. 

frame: time interval equal to 20 ms, or 160 samples at an 8 kHz sampling rate. 

harmonic noise weighting filter: this filter exploits the noise masking properties of the spectral peaks which occur at 
harmonics of the pitch frequency by weighting the residual error less in regions near the pitch harmonics and more in 
regions away from them. Note that this filter is only used when the long term filter is enabled (MODE = 1, 2 or 3). 

high pass filter: this filter is used to de-emphasize the low frequency components of the input speech signal. 

integer lags: set of lag values having whole sample resolution. 

interpolating filter: FIR filter used to estimate sub-sample resolution samples, given an input sampled with integer 
sample resolution. 

lag: long term filter delay. This is typically the pitch period, or a multiple or sub-multiple of it. 

long term filter: this filter is used to generate the periodic component in the excitation for the current subframe. This 
filter is only enabled for MODE = 1, 2 or 3. 

LPC coefficients: Linear Predictive Coding (LPC) coefficients is a generic descriptive term for describing the short 
term filter coefficients. 

open loop lag search: process of estimating the near optimal lag directly from the weighted speech input. This is done 
to narrow the range of lag values over which the closed loop lag search shall be performed. 

open loop lag trajectory: for a given frame, the sequence of near optimal lag values whose elements correspond to the 
four subframes as determined by the open loop lag search. 

reflection coefficients: alternative representation of the information contained in the short term filter parameters. 

residual: output signal resulting from an inverse filtering operation. 

short term filter: this filter introduces, into the excitation signal, short term correlation which models the impulse 
response of the vocal tract. 

soft interpolation: process wherein a decision is made for each frame to use either interpolated or uninterpolated short 
term filter parameters for the four subframes in that frame. 

soft interpolation bit: one bit code indicating whether or not interpolation of the short term parameters is to be used in 
the current frame. 

spectral noise weighting filter: this filter exploits the noise masking properties of the formants (vocal tract resonances) 
by weighting the residual error less in regions near the formant frequencies and more in regions away from them. 

subframe: time interval equal to 5 ms, or 40 samples at an 8 kHz sampling rate. 

vector quantization: method of grouping several parameters into a vector and quantizing them simultaneously. 

GSP0 vector quantizer: process of vector quantization, its intermediate parameters (GS and P0) for the coding of the 
excitation gains p and y. 
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VSELP codebook: Vector-Sum Excited Linear Predictive ( VSELP) codebook, used in the GSM half rate speech coder, 
wherein each codebook vector is constructed as a linear combination of the fixed basis vectors. 

zero input response: output of a filter due to all past inputs, i.e. due to the present state of the filter, given that an input 
of zeros is applied. 

zero state response: output of a filter due to the present input, given that no past inputs have been applied, i.e. given 
the state information in the filter is all zeroes. 



3.2 Symbols 



For the purposes of the present document, the following symbols apply: 



A(z) 
cq 

bL(n) 

P 
C(z) 

e(n) 
fjG) 

giG) 



Y 
H 



L 

^max 

^min 

M 

Ml 

M2 

MODE 

N A 



N F 

Np 

N s 

PI 

P2 

P3 

Pf 



pitch 

Ql 

Q2 

Q3 

R0 

r(n) 

rL( n ) 



Short term spectral filter. 

The LPC coefficients. 

The output of the long term filter state (adaptive codebook) for lag L. 

The long term filter coefficient. 
Second weighting filter. 
Weighted error signal 

The coefficients of the j tn phase of the 10th order interpolating filter used to evaluate candidate 
fractional lag values; i ranges from to Pf-1. 

The coefficients of the j l " phase of the 6th order interpolating filter used to interpolate C's and G's 

as well as fractional lags in the harmonic noise weighting; i ranges from to Pg-1. 

The gain applied to the vector(s) selected from the VSELP codebook(s). 

A M2 bit code indicating the vector to be selected from the second VSELP codebook (when 

operating in mode 0). 

A M or Ml bit code indicating the vector to be selected from one of the two first VSELP 

codebooks. 

The long term filter lag value. 

142 (samples), the maximum possible value for the long term filter lag. 

21 (samples), the minimum possible value for the long term filter lag. 

9, the number of basis vectors, and the number of bits in a codeword, for the VSELP codebook 
used in modes 1, 2, and 3. 

7, the number of basis vectors, and the number of bits in a codeword, for the first VSELP 
codebook used in mode 0. 

7, the number of basis vectors, and the number of bits in a codeword, for the second VSELP 
codebook used in mode 0. 

A two bit code indicating the mode for the current frame (see annex A). 

170, the length of the analysis window. This is the number of high pass filtered speech samples 

used to compute the short term filter parameters for each frame. 

160, the number of samples per frame (at a sampling rate of 8 kHz). 

10, the short term filter order. 

40, the number of samples per subframe (at a sampling rate of 8 kHz). 

6, the number of bits in the prequantizer for the rl - r3 vector quantizer. 

5, the number of bits in the prequantizer for the r4 - r6 vector quantizer. 

4, the number of bits in the prequantizer for the r7 - rlO vector quantizer. 

The order of one phase of an interpolating filter used to evaluate candidate fractional lag values. Pf 

equals 10 for j ^ and equal to 1 for j = 0. 

The order of one phase of an interpolating filter, fi(n), used to interpolate C's and G's as well as 

fractional lags in the harmonic noise weighting, P„ equals 6. 

The time duration between the glottal pulses which result when the vocal chords vibrate during 

speech production. 

1 1, the number of bits in the rl - r3 reflection coefficient vector quantizer. 
9, the number of bits in the r4 - r6 reflection coefficient vector quantizer. 

8, the number of bits in the r7 - rlO reflection coefficient vector quantizer. 
A five bit code used to indicate the energy level in the current frame. 
The long term filter state (the history of the excitation signal); n < 

The long term filter state with the adaptive codebook output for lag L appended. 
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s'(n) 
W(z) 

^-hnw 

M 
LxJ 
K 

Tx(i) 
i = j 

K 



n 



max(x,y) 
min(x,y) 
round(x) 



Synthesized speech. 

Spectral weighting filter. 

The harmonic noise weighting filter coefficient. 

The adaptive pitch prefilter coefficient. 

Ceiling function: the largest integer y where y < x + 1,0. 

Floor function: the largest integer y where y < x. 

Summation: x(j)+x(j+l)+...+x(K). 



Product: x(j)(x(j+l))....(x(K)) 

Find the larger of two numbers x and y. 

Find the smaller of two numbers x and y. 

Round the non-integer x to the closest integer y: y = I X + 0,5 ly: y=x+0,5. 



3.3 



Abbreviations 



For the purposes of the present document, the following abbreviations apply: 

AFLAT Autocorrelation Fixed point LAttice Technique 

CELP Code Excited Linear Prediction 

FLAT Fixed Point Lattice Technique 

LTP Long Term Predictor 

SST Spectral Smoothing Technique 

VSELP Vector-Sum Excited Linear Prediction 



4 Functional description of the GSM half rate speech 

codec 

The GSM half rate codec uses the VSELP (Vector-Sum Excited Linear Prediction) algorithm. The VSELP algorithm is 
an analysis-by-synthesis coding technique and belongs to the class of speech coding algorithms known as CELP (Code 
Excited Linear Prediction). 

The GSM half rate codec's encoding process is performed on a 20 ms speech frame at a time. A speech frame of the 
sampled speech waveform is read and based on the current waveform and the past history of the waveform, the codec 
encoder derives 18 parameters that describe it. The parameters extracted are grouped into the following three general 

classes: 

- energy parameters (R0 and GSP0); 

- spectral parameters (LPC and INT_LPC); 

- excitation parameters (LAG and CODE). 

These parameters are quantized into 112 bits for transmission as described in annex A and their order of occurrence 
over Abis is given in annex B. 

The GSM half rate codec is an analysis-by-synthesis codec, therefore the speech decoder is primarily a subset of the 
speech encoder. The quantized parameters are decoded and a synthetic excitation is generated using the energy and 
excitation parameters. The synthetic excitation is then filtered to provide the spectral information resulting in the 
generation of the synthesized speech (see figure 1). 
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GSM half rate speech codec 
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Figure 1 : Block diagram of the GSM half rate speech codec 

The ANSI-C code that describes the GSM half rate speech codec is given in GSM 06.06 [2] and the test sequences in 
GSM 06.07 [3] (see clause 5 for the codec homing test sequences). 

4.1 GSM half rate speech encoder 

The GSM half rate speech encoder uses an analysis by synthesis approach to determine the code to use to represent the 
excitation for each subframe. The codebook search procedure consists of trying each codevector as a possible excitation 
for the Code Excited Linear Predictive (CELP) synthesizer. The synthesized speech s'(n) is compared against the input 
speech and a difference signal is generated. This difference signal is then filtered by a spectral weighting filter, W(z), 
(and possibly a second weighting filter, C(z)) to generate a weighted error signal, e(n). The power in e(n) is computed. 
The codevector which generates the minimum weighted error power is chosen as the codevector for that subframe. The 
spectral weighting filter serves to weight the error spectrum based on perceptual considerations. This weighting filter is 
a function of the speech spectrum and can be expressed in terms of the a parameters of the short term (spectral) filter. 



W(z) 



i=i 



(i) 



The computation of the Ct x coefficients is described in subclause 4.1.7. 

The second weighting filter C(z), if used, is a harmonic weighting filter and is used to control the amount of error in the 
harmonics of the speech signal. If the weighting filter(s) are moved to both input paths to the subtracter, an equivalent 
configuration is obtained as shown in figure 2. 
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Figure 2: Block diagram of the GSM half rate speech encoder (MODE = 1,2 and 3) 

Here H(z) is the combination of A(z), the short term (spectral) filter, and W(z), the spectral weighting filter. These 
filters are combined since the denominator of A(z) is cancelled by the numerator of W(z). 



H(z) 



1 



i-E^"' 



(2) 



There are two approaches that can be used for calculating the gain, y. The gain can be determined prior to codebook 
search based on residual energy. This gain would then be fixed for the codebook search. Another approach is to 
optimize the gain for each codevector during the codebook search. The codevector which yields the minimum weighted 
error would be chosen and its corresponding optimal gain would be used for y. The latter approach generally yields 
better results since the gain is optimized for each codevector. This approach also implies that the gain term needs to be 
updated at the subframe rate. The optimal code and gain for this technique can be computed as follows. 

The input speech is first filtered by a high pass filter as described in subclause 4.1.1. The short term filter parameters are 
computed from the filtered input speech once per frame. A fast fixed point covariance lattice technique is used. 
Subclauses 4.1.3 and 4. 1 .4 describes in detail how the short term parameters are determined and quantized. An overall 
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frame energy is also computed and coded once per frame. Once per frame, one of the four voicing modes is selected. If 
MODE^O, the long term predictor is used and the long term predictor lag, L, is updated at the subframe rate. L and a 
VSELP codeword are selected sequentially. Each is chosen to minimize the weighted mean square error. The long-term 
filter coefficient, p, and the codebook gain, y, are optimized jointly. Subclause 4.1.8 describes the technique for 
selecting from among the voicing modes and, if one of voiced modes is chosen, determining the long-term filter lag. 
Subclause 4.1.10 describes an efficient technique for jointly optimizing p, y and the codeword selection. 
Subclause 4.1.10 also includes the description of the fast VSELP codebook search technique. The p and y parameters 
are transformed to equivalent parameters using the frame energy term, and are vector quantized every subframe. The 
coding of the frame energy and the p and y parameters is described in subclause 4.1.11. 

4.1.1 High-pass filter 

The 13 bit linear Pulse Code Modulated (PCM) input speech, x(n), is filtered by a fourth order pole-zero high pass 
filter. This filter suppresses the frequency components of the input speech which are below 120 Hz. The filter is 
implemented as a cascade of two second-order Infinite Impulse Response (IIR) filters. Incorporated into the filter 
coefficients is a gain of 0,5. The difference equation for the first filter is: 

2 2 

y(n) = y^jxin - i) + 2_, a l, jy\ n ~ J) (3) 

(=0 j=\ 

where: 

b 10 = 0,335052 

b 1 1 = -0,669983 a 1 1 = 0,926 1 1 7 

bi2 = 0,335052 ai2 = -0,429413 

The difference equation for the second filter is: 

2 2 

y(n) = ^b 2i y(n-i) + Y, a 2jy{ n ~ j) ( 4 ) 

(=0 ;=1 

where: 

b20 = 0,335052 

b 2 i = -0,669434 a21= 0,965332 

b 2 2 = 0,335052 a22 = -0,469513 



4.1.2 Segmentation 

A sample buffer containing the previous 195 input high pass filtered speech samples, y(n), is shifted so that the oldest 
160 samples are shifted out while the next 160 input samples are shifted in. The oldest 160 samples in the buffer 
correspond to the next frame of samples to be encoded. The analysis interval comprises the most recent 170 samples in 
the buffer. The samples in the buffer are labelled as s(n) where < n < 194 and s(0) is the first (oldest) sample. 

4.1 .3 Fixed Point Lattice Technique (FLAT) 

Let r; represent the j tn reflection coefficient. The FLAT algorithm for the determination of the reflection coefficients is 
stated as follows: 
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STEP 

1 



Compute the covariance (autocorrelation) matrix from the input speech: 



STEP 

2 



STEP 

3 



STEP 

4 

STEP 

5 



N A 

(j{i,k)= ^ sin + 24 - i)s(n + 24 - k) < i, k< 

n=N p 

The (|)(i,k) array is modified by windowing 
0'(i,k) = 0(i,k)w(\i-k\) 

F (i,k) = </>'(i,k) 
B (i,k) = f{i + l,k + l) 
C o (i,k) = 0'(i,k + l) 

setj = 1 



N r 



(5) 



< i, k < N p 


(6) 


< i, k < N p - 1 


(7) 


< i, k < N p - 1 


(8) 


< i, k < N p - 1 


(9) 



Compute rj 



Cj.M + Cj^Np-j^p 



STEP 
6 



^_l(0,0) + Bj.^Ofi) + Fj^Np - j,N p - j) + Bj.^Np - j, N p - j 
If j = Np then done. 



(10) 



STEP 

7 



Update Fj(i,k), Bj(i,k), Cj(i,k) < i, k < N P -j-l 

Fj{i,k) = Fj^ii^ + r^Cj^iu^ + Cj^ikJjj + rJBj^iuk) (11) 

Bj{i,k) = B hl {i + \,k + l) + r^Cj.^i + l,k + l) + C hl {k + \,i + l)) + rJFj^i + l,k + l) , (12) 

Cj {i,k)= Cj_ x {i, k + l) + ^ [b hl {i, k + l) + Fj_ { {i,k + l)) + rj C hl {k + 1, i) ( 1 3) 

STEP j=j+l 



STEP go to step 5. 
9 

The windowing coefficients, w(li-kl), are found in the table 1. 

Table 1 : Windowing coefficients 



w(0) 


0,998966 




w(5) 


0,974915 


w(1) 


0,996037 


w(6) 


0,969054 


w(2) 


0,991663 


w(7) 


0,963060 


w(3) 


0,986399 


w(8) 


0,956796 


w(4) 


0.980722 


w(9) 


0,950127 
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This algorithm can be simplified by noting that the 0', F and B matrices are symmetric such that only the upper 
triangular part of the matrices need to be computed or updated. Also, step 7 is done so that Fj(i,k), Bj(i-l,k-l), Cj(i,k-1), 
and Ci(k,i-1) are updated together and common terms are computed once and the recursion is done in place. 

4.1.4 Spectral quantization 

A three segment vector quantizer of the reflection coefficients is employed. A reduced complexity search technique is 
used to select the vector of reflection coefficients for each segment. The reflection coefficient vector quantizer 
codebooks are stored in compressed form to minimize their memory requirements. 

The three segments of the vector quantizer span reflection coefficients r j ri>, r4 i^, and ij - r j q respectively. The bit 
allocations for the vector quantizer segments are: 

Q4 11 bits 

Q 2 9 bits 

Q 3 8 bits 

A reflection coefficient vector prequantizer is used at each segment. The prequantizer size at each segment is: 

Pi 6 bits 

P 2 5 bits 

P 3 4 bits 

At a given segment, the residual error due to each vector from the prequantizer is computed and stored in temporary 
memory. This list is searched to identify the four prequantizer vectors which have the lowest distortion. The index of 
each selected prequantizer vector is used to calculate an offset into the vector quantizer table at which the contiguous 
subset of quantizer vectors associated with that prequantizer vector begins. The size of each vector quantizer subset at 
the k-th segment is given by: 

2 Qk 

S k =-]T- (14) 

2 k 

The four subsets of quantizer vectors, associated with the selected prequantizer vectors, are searched for the quantizer 
vector which yields the lowest residual error. Thus at the first segment, 64 prequantizer vectors and 128 quantizer 
vectors are evaluated, 32 prequantizer vectors and 64 quantizer vectors are evaluated at the second segment, and 16 
prequantizer vectors and 64 quantizer vectors are evaluated at the third segment. 

4.1 .4.1 Autocorrelation Fixed Point Lattice Technique (AFLAT) 

An autocorrelation version of the FLAT algorithm, AFLAT, is used to compute the residual error energy for a reflection 
coefficient vector being evaluated. Compute the autocorrelation sequence R(i), from the optimal reflection coefficients, 
r;, over the range < i < Np. 

STEP Define the initial conditions for the AFLAT recursion: 
1 

P Q {i) = R{i), 0<i<N p -i (15) 

V Q {i) = R(\i + l\), l-N p <i<N p -l (16) 

STEP Initialize k, the vector quantizer segment index: 
2 

k=l (17) 

STEP Let Ij(k) be the index of the first lattice stage in the k-th segment, and I n (k) be the index of the last 

3 lattice stage in the k-th segment. 
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STEP Initialize j, the index of the lattice stage, to point to the beginning of the k-th segment: 

4 

J = lM (18) 

STEP Set the initial conditions Pj.j and Vj_i to: 

5 

P j . l {i) = P hl {i), 0<i<I h {k)-I,{k) (19) 

V j _ l {i) = V H {i), -7 ft (fe) + / / (fe)</</ A (fe)-/ / (fe) (20) 

STEP Compute the values of Vj and Pi arrays using: 
6 



Pj (i) = (l + ff)Pj_ x (i) + fj [F 7 _! (j) + V M (- i)\, 0<i<I h (k)-j-l (21) 

Vj (i) = Vj_i (i + 1) + ffv H (- i - 1) + 2r,- P^j ( \i + 1|), l + j-N p <i<N p -j-l (22) 

STEP Increment j: 

7 

j=j+l 

STEP If j < Ih(k) go to STEP 6. 
8 

STEP The residual error out of lattice stage Ih(k), given the reflection coefficient vector r , is computed using 
" equation (21): 

E r = P h{k] (0) (23) 

STEP Using the AFLAT recursion outlined, the residual error due to each vector from the prequantizer at the 

10 k-th segment is evaluated, the four subsets of quantizer vectors to be searched are identified, and 
residual error due to each quantizer vector from the selected four subsets is computed. The index of r , 
the quantizer vector which minimized E r over all the quantizer vectors in the four subsets, is encoded 
with Qk bits. 

STEP If k < 3 then the initial conditions for doing the recursion at segment k+1 need to be computed. Set j, 

1 1 the lattice stage index, equal to: 

j = ll{k) (24) 

STEP Compute: 
12 

Pj (j) = (l + rf )Pj_ x (i) + 7j Vj_ x (i) + Vj_ x (- i) , 0<i<N p -j-l (25) 

Vj (i) = Vj_i (i + 1) + rfVj_i (- i - 1) + If ) Pj_ r {\i + 1|), 1 + j - N p < i < N p - j - 1 (26) 



STEP Increment j , 
13 

j=j+l 

STEP If j < I h (k) go to STEP 12 

14 

STEP Increment k, the vector quantizer segment index: 

15 

k=k+l 
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STEP Ifk<3gotoSTEP4. 
16 

Otherwise, the indices of the reflection coefficient vectors for the three segments have been chosen, and 
the search of the reflection coefficient vector quantizer is terminated. 

To minimize the storage requirements for the reflection coefficient vector quantizer, eight bit codes for the individual 
reflection coefficients are stored in the vector quantizer table, instead of the actual reflection coefficient values. The 
codes are used to look up the values of the reflection coefficients from a scalar quantization table with 256 entries. 

4.1 .5 Frame energy calculation and quantization 

The unquantized value of R0, R(0), is computed during the computation of the short term predictor parameters. 

^(0.0)^(10.10) 

v ' 320 

where (|)(i,k) is defined by equation (5). R(0) is then converted into dB relative to full scale (full scale, Rmax> is defined 
as the square of the maximum sample amplitude). 



^B=101og 10 



( R[0) ^ 



V "max J 



(28) 



RdB is then quantized to 32 levels. The 32 quantized values for R^B range from a minimum of -66 (corresponding to a 
code of for R0) to a maximum of -4 (corresponding to a code of 3 1 for R0). The step size of the quantizer is 2 (2 dB 
steps). R0 is chosen as: 

R0 which minimizes abs(R0 - (RdB + 66)/2) (29) 

where R0 can take on the integer values from to 3 1 corresponding to the 32 codes for R0. 

Decoding of the R0 code is given by: 

R(o) = R m ,M <!RO) - 66)m m 

4.1 .6 Soft interpolation of the spectral parameters 

Interpolation of the short term filter parameters improves the performance of the GSM half rate encoder. The direct 
form filter coefficients (OCj's), which correspond to quantized reflection coefficients, are the spectral parameters used for 
interpolation. The GSM half rate speech encoder uses either an interpolated set of OCj's or an uninterpolated set of OCj's, 
choosing the set which gives better prediction gain for the frame. 

Two sets of LPC coefficient vectors are generated: the first corresponds to the interpolated coefficients, the second to 
the uninterpolated coefficients. The frame's speech samples are inverse filtered using each of the two coefficient sets, 
and the residual frame energy corresponding to each set is computed. The coefficient set yielding the lower frame 
residual energy is then selected to be used. If the residual energies are equal, the uninterpolated coefficient set is used. 
INT_LPC, a soft interpolation bit, is set to 1 when interpolation is selected or to otherwise. 

To generate the interpolated coefficient set, the coder interpolates the OCj's for the first, second, and third subframes of 
each frame. The fourth subframe uses the uninterpolated OCj's for that frame. 

The interpolation is done as follows. Let aj l be the direct-form LPC coefficients corresponding to the last frame, (X[ q 
be the direct-form LPC coefficients corresponding to the current frame, and Del to be the interpolation curve used. The 
interpolated direct-form LPC coefficient vector at the j-th subframe of the current frame, 0Cj ;, is given by: 

a i j=a iL +Del(j,INT_SOFT)(a ic -a iL ), l < i < Np, l < j < 4 (31) 

The values of the interpolation curve Del are given in table 2. 
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Table 2: Values of the interpolation curve Del 



j 


Del(j,0) 


Del(j,1) 


1 


0,0 


0,30 


2 


1,0 


0,62 


3 


1,0 


0,92 


4 


1,0 


1,00 



From this point on, the subframe index j is omitted for simplicity when referring to an coefficients, although it is 
implied. For interpolated subframes, the oci's are converted to reflection coefficients to check for filter stability. If the 
resulting filter is unstable, then uninterpolated coefficients are used for that subframe. The uninterpolated coefficients 
used for subframe 1 are the previous frame's coefficients. The uninterpolated coefficients used for subframes 2, 3, and 4 
are the current frame's coefficients. 



4.1 .7 Spectral noise weighting filter coefficients 

To exploit the noise masking potential of the formants, spectral noise weighting is applied. The computation of the CC; 

coefficients, used by spectral noise weighting filters W(z) and H(z), is now described. Define an impulse sequence 8(n) 
over N s samples: 



(5(0) = 1,0 
S(n) = 0,0 



(32) 



where 1 < n < N s -1 and li3(n) is the zero-state response of the cascade of three filters to 8(n). The three filters are an 
LPC synthesis filter, an inverse filter using a weighting factor of 0,93 and a synthesis filter with a weighting factor of 
0,7. In equation form: 



Np 



h { (n) = S(ri) + ^ a t \ (n-i) 



i=\ 



Np 



h 2 (n) = h i (n)-Y J (0,93) i a i h l (n-i) 



i=\ 



N P 



h 3 (n) = h 2 (n) + Y J (0J) i a i h 3 (n-i), 



0<n<N s -l 



0<n<N s -l 



0<n<N s -l 



(33) 



(34) 



(35) 



/=1 



where oq's are the direct form LP coefficients. The autocorrelation sequence of h3(n) is calculated using: 



Ns-1 

n=i 



< i < N r 



(36) 



From Rfo (i) the reflection coefficients which define the combined spectrally noise weighted synthesis filter are 
computed using the AFLAT recursion once per frame. 



STEP 1 Define the initial conditions for the AFLAT recursion: 

P (i) = R h3 (i), 

V (i) = R h3 {\i + l\), 

STEP 2 Initialize j, the index of the lattice stage, to point to the first lattice stage: 

j = l 



< i < N p -1 



l-N p <i<N p -l 



(37) 



(38) 
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STEP 3 



Compute r;, the j-th reflection coefficient, using: 

vj-m 



r j = 



0-i(°) 



(39) 



STEP 4 



Given r;, update the values of Vi and Pi arrays using 



J 



r J 



0<i<N p -j-l (40) 
l+j-N p <i<N p -j-l(41) 



STEP 5 

STEP 6 

STEP 7 



Pj (i) = {l + rf) Pj_ { (i) + rj [Vj_! (/) + Vj_! (-/) 

Vj (/) = V H (i + 1) + rj V M (-/ - 1) + 2rj Pj_ x ( \i + 1| 
Increment j : 

If j < Np go to STEP 3, otherwise all Np reflection coefficients have been obtained. 

The reflection coefficients, rj, are then converted to direct-form LPC filter coefficients,^ for 
implementing the combined spectrally noise weighted synthesis filter H(z) and the filter W(z). 



The method for the spectral noise weighting filter coefficient update mimicks how the direct form LPC filter 
coefficients are updated at subframes of a frame (subclause 4. 1 .6). No stability check of interpolated spectral noise 
weighting filter coefficients is done at subframes 1, 2, or 3 if the interpolation flag, INT_LPC="1", but if uninterpolated 
coefficients are used at subframes 1, 2, and/or 3 due to instability of the unweighted coefficients (INT_LPC = "0"), 
uninterpolated weighting filter coefficients are also used at those subframes. 

4.1 .8 Long Term Predictor lag determination 

Figure 3 illustrates that the long term lag optimization looks just like a codebook search where the codebook is defined 
by the long term filter state and the specific vector in the codebook is pointed to by the long term predictor lag, L. The 
input p(n) is the weighted input speech for the subframe minus the zero input response of just the H(z) filter. 



P 



X b(n 



P(n) 



Long term 
filter state 



A- 



H(z) 



^-k^ 1 * 



I 



total weighted 
error 



Figure 3: Long term predictor lag search 



The GSM half rate speech encoder uses a combination of open loop and closed loop techniques in choosing the long 
term predictor lag. First an open loop search is conducted to determine "candidate" lags at each subframe. Then at most, 
two best candidate lags at each subframe are selected, with each serving as an anchor point for constructing an open 
loop frame lag trajectory, subject to a maximum delta coding constraint. The frame lag trajectory which minimizes the 
open loop LTP spectrally weighted error energy for the frame is then chosen. The open loop LTP prediction gains 
corresponding to the winning trajectory are used to select the voicing mode 1, 2 or 3. If MODE^O, the closed loop lag 
evaluation is initiated. The winning trajectory has associated with it a list of lags to be searched closed loop at each 
subframe. 

It is possible to allow L to take on fractional values, thus increasing the resolution, and in turn the performance, of the 
adaptive codebook. Table 3 shows the allowable lags. 
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Table 3: Allowable lags 



Range 


Resolution 


Number of lags in range 


21 to 22 2/3 


1/3 


6 


23 to 34 5/6 


1/6 


72 


35 to 49 2/3 


1/3 


45 


50 to 89 1/2 


1/2 


80 


90 to 142 


1 


53 



The resolution of the long term filter state may be increased by upsampling and filtering the state. In this 
implementation, a non-causal, zero-phase Finite Impulse Response (FIR) filter is used. Where needed, the future 
samples for the non-causal filtering operation are replaced by the output of the predictor. 



4.1.8.1 



Open loop long term search initialization 



An open-loop lag search is done to narrow the range of lags over which a closed-loop search will eventually be 
performed. 

The first steps of the open-loop subframe lag search are as follows: 

STEP 1 Initialize the subframe counter 

m=l 

STEP 2 The autocorrelation sequence of y(n) the input speech, s(n), filtered by W, is calculated for all allowable 
integer lags, and for a few integer lags below and above the lower and upper limits for the current 
subframe. 

tf,-l p p 

C(k,m)= ^y(« + (m-l)/vJy(n + (m-l)/V,-4 L^ n - -±- < k < L max + -^-- 1 (42) 

where L m j n = 21 and L max =142. 

The value P„ is the order of one phase of the interpolating FIR filter used to interpolate the correlations. 
The energy of y(n) for the subframe is computed: 

N s -l p p 

G(k,m)=Y,y 2 (n + {m-l)N s -k), L Mn -- 1-< k < L max + y~l (43) 
n=0 

S 1 EP 3 These arrays, C(k,m) and G(k,m), are searched for the integer lag which maximizes C^(k,m)/G(k,m) 
where C(k,m) and G(k,m) need to be greater than 0. 

STEP 4 If a valid maximum is found in step 3, the values for the lag, C, and G corresponding to the maximum are 
retained in the arrays as Lp ea k(0,m), Cp ea k(0,m), and Gp ea k(0,m). 



Otherwise, L pea k(0,m) = L m i n 
C peak (0,m) = 

G peak(°> m ) = 1 
STEP 5 m=m+l 
STEP 6 If m < 4, go to step 2 



(44) 



(45) 
(46) 
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STEP 7 Calculate the open loop frame LTP prediction gain: 

4 



P v = 10 log 



10 



£*(0,m) 



m=\ 



m=\ 



R(0,m)- 



G peak (0,m) 



(47) 



where 



N.-l 



R(0,m)= ^y 2 (n + {m-l)N s ). 



n=0 



(48) 



STEP 8 Determine if the voicing mode is unvoiced: 

If P v < 1,7 then MODE=0, the long term predictor is disabled and the open loop and closed loop lag 
searches are aborted. In this case, proceed to subclause 4.1.10. 



4.1.8.2 



Open loop lag search 



When MODE ■£ 0, the lag search processing is continued. The next part of the search finds the allowable lag (see 

c 2 

table 3) which maximizes in the vicinity of the best open-loop integer resolution lag, Lp ea k(0,m), for values of 

G 
C>0. 

STEP 1 Initialize the subframe counter 

m=l 

Initialize the peak index 



STEP 2 



STEP 3 



^,m _ U 

Using interpolated versions of the C and G arrays, allowable lag values k' in the range: 

Lpeak(0,m) - 1 < k < L peak (0,m) + 1 (49) 

are searched for a k which maximizes 



Cf(k) 
Gt(k) 

where 



(50) 
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STEP 4 



STEP 5 



STEP 6 



STEP 7 



C I (k) = Y J 8j(i)C{[^-3 + i,m) 
i=0 

5 
G I (k) = Y J gjd)G{[^-3 + i,m) 



(51) 



(52) 



=0 



and 

j = 6{\k-]-k) 



(53) 



The coefficients of the interpolating filter are g;(i) for < i < 5. 

Only Q(k) > and Gi(k) > values are considered. If no positive correlation is found, then set 

^hnw,m = 0, L peak (l,m)=L min , and go to Step 22. 

Otherwise, store the information related to the valid best allowable lag k. 

(54) 



Ep,m — Ep,m + 1 



Lpeak(Lp,m> m )-k 

Cpeak(Lp,m,m)=Ci(k) 

Gpeak(Lp,m,m)=Gi(k) 



(55) 



(56) 
(57) 



C l 



The next part of the search evaluates , for C > and G > 0, at the submultiples of the lag 

G 
Lpeak(Lp,m>m) to find candidate peaks. 

Initialize the divisor 

J = 2 

Find nearest integer lag corresponding to submultiple of maximum peak 

k 1 =round[L p eak(l,m)/J] (58) 

Determine if submultiple is within allowable lag range 

If ki < Lmin 

Go to step 12 

Find value of k where C^(k',m)/G(k',m) is a maximum for 

max(L m i n ,ki-3) < k< min(L max ,ki+3) (59) 

If either C(k',m) < or G(k',m) < go to step 11. 
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STEP 8 



STEP 9 



STEP 10 



Determine if maximum in step 7 is a peak 


If 




C 2 (k'-l,m) ^ 
G(k'-l,m) 


C (k',m) 

G(k' ,m) 


Go to step 1 1 




If 




C 2 {k'+l,m) 

> 

G{k'+\,m) 


C \k\m) 
G{k' ,m) 


Go to step 1 1 





(60) 



(61) 



A peak has been found at an integer lag, k. Using interpolated versions of the C and G arrays, 
allowable lag values within + 1 (exclusive) of k are searched. 

Find k where 



Gj(k) 

is a maximum, where 

5 

C I (k) = Yjgj<J)c{\k~\-3 + i,m) 

5 
G I {k) = Y J gj(i)G{\k\-3 + i,m) 

where 

j = 6(|Y|-*) 

and k'-l <k<k'+l 

Only Cj(k) > and Gj(k) > are considered. 



(62) 



(63) 



(64) 



(65) 



(66) 



If the prediction gain exceeds a threshold, the corresponding lag, C\, and G\ are stored in the 
Lpeak(X CpeakO, and Gp ea k() arrays; otherwise, these values are not stored. 

If 



Cf(k) 

G r (k) 



>R(0,m) 



R(0,m) 
10 x 



where x = 7,5 log 



10 



R{0,m) 



R{0,m)- 



C 



peak 
'peak 



(0,m) 
(0,m) 



(67) 
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STEP 1 1 



STEP 12 



STEP 13 



STEP 14 



then 



Ep,m— Ep,m+1 



(68) 

(69) 
(70) 
(71) 



Lpeak(Lp,m> m )-k 
C p eak(L p ,m,m)=Ci(k) 

Gpeak(L p ,m,m)=G I (k) 

Increment divisor and check the next submultiple 

J=J+1 

Go to step 5 

A full-resolution search (1/6 sample resolution) is done for a peak within 1 integer lag (exclusive) 
of the shortest lag. 

Find k such that 



Cj(k) 
Gt(k) 

is a maximum, where 

5 



(72) 



C l (k) = Y J gjd)C{[^-3 + i,m) 



;=0 

5 



G I {k) = Y J gj(i)G{\k\-3 + i,m) 



(73) 



(74) 



(=0 



(75) 



j = b{\k-]-k) 

L xmn-^ L peak^ L p,m)- 1 ) < k < min [ L max + g. L peak { L p,m) + l j ( 7 6) 

The fractional lag corresponding to the maximum is referred to as Lpjtch,m- This lag is used by the 
harmonic noise weighting function C(z) at subframe m. Then 

^ pitch. m ~ *- I\ pitch. mj ( ' ' ) 

" pitch, m = Gj y ^pitch,m J ( ' °) 

The harmonic noise weighting coefficient for subframe m is calculated in this step (see 
subclause 4.1.9) 



c, 






u pitch, m 



(79) 



Once all the correlation peaks associated with sub multiples of the Lp ea k(l,m) have been examined, 
the correlation peaks associated with multiples of Lpit cn ^ m are examined. 



Initialize the multiplier 
J = 2 
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STEP 15 



STEP 16 



STEP 17 



STEP 18 



STEP 19 



Find nearest integer lag corresponding to a multiple of the fundamental lag 

kj = round [Lpi tc h,m* J ] ( 80 ) 

Determine if multiple is within allowable lag range 

If 

M> E max 

Go to step 22 

Find value of k' where C^(k',m)/G(k',m) is a maximum for 

max(L m i n ,ki-3) < k '< min(L max ,ki+3) 

If either C(k',m) < or G(k',m) < go to step 21. 

Determine if maximum in step 17 is a peak 

If 



(81) 



C z (fc'-l,m) C z (k\m) 
G[k'~ l,m) G(k\m) 



Go to step 21 
If 



C 2 {k'+\,m) C 2 (k\m) 



(82) 



G[k'+\,m) G(k',m) 
Go to step 21 



(83) 



A peak has been found at an integer lag, k. Using interpolated versions of the C and G arrays, 
allowable lag values within + 1 (exclusive) of k are searched. 

Find k where 



is a maximum, where 

5 



(84) 



C I (k) = Y J gj(i)C{\k~]-3 + i,m) 
i=0 

5 
G I (k) = Y J gjd)G{[^-3 + i,m) 
i=0 

where 

./ = 6flY|-*) 

and k'-l <k<k'+l 

Only Ci(k) > and Gi(k) > are considered. 



(85) 



(86) 



(87) 



(88) 
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STEP 20 If the prediction gain exceeds a threshold, the corresponding lag, Cj, and G\ are stored. 

If 



CJ(k) _.. , i?(0,m) 
' > R(0,m) , where x = 7,5 log 10 



Gj(k) 



1(T 



fl(o,m) 



fl(o,m)- 



C 



peak 



(0,m) 



1 peak 



(0,m) 



(89) 



then 



STEP 21 



STEP 22 



4.1.8.3 



Ep,m — Ep,m + I 










(90) 


Epeak(Ep,m) = k 










(91) 


Cpeak(L p ,m)=Cl(k) 










(92) 


Gpeak(L p ,m)=Gl(k) 










(93) 


Increment multiplier 


and check the nex 


: multiple 




J=J+1 












Go to step 15 












Increment subframe 


pointer 


and 


repeat 


for all subframes 


m=m+l 












If m<4 












Go to step 2. 













Otherwise, the list of correlation peaks and the harmonic noise weighting filter parameters for each 
subframe have been found. 



Frame lag trajectory search (Mode * 0) 



The frame lag trajectory search uses the list of potential lag values to determine the one lag value for each subframe 
which minimizes the open loop prediction error energy for the frame subject to the constraints of the delta lag coding 
employed for subframes 2, 3 and 4. Several candidate lag trajectories are determined. The trajectory which minimizes 
the open loop prediction error energy for the frame is chosen. 

c 2 

In subclause 4.1.8.2, the open loop lag search found a list of lags, Lp ea k(i,m), corresponding to the peaks, for each 

G 

C 2 

subframe. Each trajectory evaluation begins with one of the subframes and selects a lag corresponding to a peak 

G 
for that subframe as the anchor for that candidate trajectory. 

A maximum of 2 trajectories are anchored per subframe. From the anchor lag, the trajectory is extended forward and 
backward to the adjacent subframes in the frame subject to the lag differential coding constraints. The lag for each 
subframe on the trajectory is chosen to minimize the open loop frame prediction error energy. The trajectory search is 
described below. 

The steps involved in the frame lag trajectory evaluation and selection are: 

STEP 1 Set m, the pointer to the selected subframe, equal to 1. 
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STEP 2 Choose the lag at the selected subframe, m, to be an anchor lag for the frame lag trajectory; i.e., the 
frame lag trajectory being evaluated needs to pass through that lag. The lag which is chosen, 

Q 2 

corresponds to the highest peak in the list of peaks at subframe m, which has not been crossed 

Gj 

by a trajectory evaluated previously. If no peaks qualify, no peaks are left, or two trajectories have 
already been anchored and evaluated at subframe m, go to step 7. Otherwise, compute the open loop 
subframe weighted error energy corresponding to the chosen lag, and store the result in the frame 
weighted error accumulator corresponding to the trajectory currently being evaluated. 

STEP 3 If m < 4, begin the forward search: 

STEP 3a Define the current subframe to be m+1. 

STEP 3b Define the forward search range as -7 to +6 levels relative to the current 

subframe's lag level. 

STEP 3c Check that the lower bound does not point to a level below the lowest allowable 

lag level, clipping if necessary. Similarly, check that the upper bound does not 
point past the highest allowable lag level; clip if necessary. 

STEP 3d Cj 

Find the lag within the range which maximizes 



NOTE: negative values of Cj are allowed. Compute the open loop subframe weighted error energy corresponding 
to that lag at the current subframe, and add the result to the frame weighted error accumulator 
corresponding to the trajectory being evaluated. 

STEP 3e If the current subframe < 4, increment the pointer to the current subframe, and go 

to step 3b. 

STEP 4 If m > 1, initiate the backward search: 

STEP 4a Define the current subframe to be m-1. 

STEP 4b Define the backward search range as -6 to +7 levels relative to the current 

subframe's lag level. 

STEP 4c Check that the lower bound does not point to a level below the lowest allowable 

lag level, clipping if necessary. Similarly, check that the upper bound does not 
point past the highest allowable lag level; clip if necessary. 

STEP 4d C, 

Find lag within the range which maximizes ' 



IGj 

NOTE: negative values of C\ are allowed. Compute the open loop subframe weighted error energy corresponding 
to that lag at the current subframe, and add the result to the frame weighted error accumulator 
corresponding to the trajectory being evaluated. 

STEP 4e If the current subframe index is > 1, decrement the pointer to the current subframe, 

and go to step 4a. 

STEP 5 Store the lags defining the frame lag trajectory derived and the open loop LTP frame weighted error 
energy which this trajectory yields. Increment the counter of evaluated frame lag trajectories. 

STEP 6 Go to step 2. 

STEP 7 If m < 4, increment m and go to step 2. 

STEP8 Choose, from the set of constructed frame lag trajectories, a lag trajectory which yields the lowest 
LTP weighted error energy for the frame, as the selected frame lag trajectory. 
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4.1.8.4 



Voicing mode selection 



The frame lag trajectory is specified by a vector K={ki,k2,k3,k4.}, where k m is the open loop LTP lag at the m-th 
subframe. Define the interpolated correlation of the input spectrally weighted speech y(n) at the m-th subframe, 
specified by lag k m , as Ci(k m ,m) and the interpolated energy of y(n), delayed by k m samples relative to the m-th 
subframe, as Gi(k m ,m). 

The open loop LTP prediction gain in dB at the m-th subframe is: 



10 log 



10 



R{0,m) 



R{0,m)- 



CJ{k m ,m) 



(94) 



Gi[k m ,m) _ 

The open loop frame LTP prediction gain, is given by: 
4 



10 log 



10 



£*(0,m) 



m=\ 



m=\ 



/?(o,m)- 



(0,m) 



' peak 



(95) 



The rules for mode selection are specified as follows: 



MODE=0ifP v < 1,7 



MODE=l if P v > 1,7 and Pm < 3,5 for any m 
MODE=2 if P m > 3,5 for all m and Pm < 7 for any m 
MODE=3 if P m > 7,0 for all m 



(96) 
(97) 
(98) 
(99) 



4.1.8.5 



Closed loop lag search 



From the selected frame lag trajectory, develop a list of lags to be searched closed loop. At each subframe, three 
allowable lag levels centered around the subframe lag, specified by the selected frame lag trajectory, will be searched. If 
the lag points to the lowest or the highest level in the table of quantized lag values, only two closed loop lag evaluations 
will be done at that subframe, with the lag outside the quantizer range being eliminated from consideration. The closed 
loop evaluation of the subframe lags is not performed if MODE=0. What follows is a description of the construction of 
the output of the long term predictor (adaptive codebook) for a given, possibly fractional lag, L. Defining: 



Lmax 

r(n) 

rL(n) 

bL(n) 

Pf 

fj(i) 

N s 



maximum possible value for long term lag L 

long term filter state; n < (history of the excitation signal) 

long term filter state with adaptive codebook output for L appended 

output of long term filter state (adaptive codebook) for lag L 

order of one phase of the interpolating FIR filter (Pf =10 except for the special case when j = 0, 
see below) 

coefficients of jth phase of interpolating FIR filter, i=0 to i=Pf -1 
number of samples per subframe (N s = 40) 



The sequence rjjn) is defined as: 



ETSI 



3GPP TS 46.020 version 7.0.0 Release 7 



28 



ETSI TS 146 020 V7.0.0 (2007-06) 



(») = 



p,-i 



"•(«) 



2 fSV L 



i = 



n-A — +/ 

2 



^max^- 1 



;0<n<N s -l 



(100) 



where; q = 



n + L + - 
6 



L, A = |_gJ, and j = 6[q-\_q\) 



The portion of the sequence rL(n) from n=0 to n=N s -l shall be calculated in order from to N s -1, so that the necessary 
terms in the sum will be available. The tn phase of the interpolating filter, f (i) , is a special case and has only one 
non-zero tap, so that if q is an integer, the summation reduces to the single term, rL(n-q). 

The output of the codebook for lag L is just the last N s samples in the sequence rL(n). 

b L (n) = r L (n) ;0<n<N s -l (101) 

C 2 

The closed loop search minimizes the weighted error by maximizing the term , where 

G 



N s -1 
c = X b L (n)p(n) 

n=0 

G= ^h (n) 

n=0 



(102) 



(103) 



The sequence b'L(n) is the zero state response of H(z) to the adaptive codebook output for lag L. The sequence p(n) is 
the input speech, weighted by the filter W(z), minus the zero input response of H(z). The error minimization is done 

c 

over only those lags in the list supplied by the open loop search. The lag L which maximizes , — (C is allowed to be 

yjG 
negative) is then chosen as the lag for the subframe. 



4.1 .9 Harmonic noise weighting 

If MODE = 1, 2 or 3, then C(z), the harmonic noise weighting transfer function, is activated. The excitation codebook 
vector and gains codebook vector are selected to minimize the spectrally and harmonically weighted error. The 
harmonic weighting filter, C(z), can be expressed as: 



-L n 



C(z) = l-A hnw z - pitch 



where 



c„ 



pilch 



(104) 



(105) 



Cpitcfb Gpitch ar, d Lpitch were determined during the open loop lag search where the subscript m denoting the 
subframe has been dropped from Lpjtch,m an d ^-hnw>m f° r notational convenience. Lpitch can ta k e on fractional values 
so the interpolating filter employed for the open loop lag search is utilized to generate the fractionally delayed samples. 
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Let x(n) represent the input of the harmonic noise weighting filter and y(n) represent the output. The filter can then be 
described by equation (106) and equation (107). 



y(n) = x(n) - A-hnw x(n - L p i tc h) 



(106) 



tin- 



P.-i 



lpitch 



2l*j®' 



lpitch 



+ 1 



(107) 



[=0 



where 



\ L pitch ~ \_ L pitch \p 



Figure 4 incorporates harmonic noise weighting (for MODE = 1, 2 or 3) and shows the VSELP excitation source. All 
error minimizations done after lag selection utilize the combination of spectral and harmonic noise weighting. 

For MODE £ 0, P(n) is the input speech signal weighted by W(z)C(z) minus the input response of H(z)C(z). For MODE 
= 0, P(n) is the input speech signal weighted by W(z) minus the zero input response of H(z). 
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Figure 4: Long term predictor lag and code search 

The zero state response of H(z)C(z) to the pitch predictor (adaptive codebook) vector, t>L(n), needs to be computed 
prior to the search of the excitation codebook. This spectrally and harmonically weighted adaptive codebook vector is 
represented by b"L(n). 



4.1 .1 Code search algorithm 

For MODE^O, the excitation codebook search procedure takes place after the long term predictor lag, L, has been 
determined. The codebook search procedure then chooses one codevector from the VSELP codebook. The GSM half 

rate speech encoder uses an excitation codebook of 2^ codevectors which is constructed from M basis vectors. 



nth 



Defining v m (n) as the m ln basis vector and uj(n) as the i m codevector in the codebook, then 



th, 



M 



i.(n)= Ye- V (n) 
m — \ 



(108) 
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where < i < 2-1; < n < N s -1. In other words, each code vector in the codebook is constructed as a linear 
combination of the M basis vectors. The linear combinations are defined by the 8 parameters. 

0i m is defined as: 

0j m = + 1 if bit m of codeword i = 1 

0j m = - 1 if bit m of codeword i = 

The codebook construction for the GSM half rate speech encoder can be restated as follows. Codevector i is constructed 
as the sum of the M basis vectors where the sign (plus or minus) of each basis vector is determined by the state of the 
corresponding bit in codeword i. The codebook search procedure finds the codevector which will produce the minimum 
total spectral and harmonic weighted error for the subframe given b"L(n) (the zero state response of H(z)C(z) to bL(n) ) 

and allowing both the gain, y, and the long term filter coefficient, (3, to be optimized for each codevector being 
evaluated. 

The filtered codevector fj(n), can be expressed as: 

M 

fiW = Ii mflm<ri (109) 

m=l 
where qm( n ) is the zero state response of H(z)C(z) to basis vector v m (n). 

If MODE=0, two VSELP codebooks are the excitation sources and are searched sequentially to identify the codeword I 
specifying the codevector selected from the first codebook, and codeword H, identifying the codevector chosen from 
the second VSELP codebook. Harmonic noise weighting is not used. When searching the second VSELP codebook, 
each codevector is evaluated assuming optimal gains for the codevector I, and the potential codevector H. 

4.1 .1 0.1 Decorrelation of filtered basis vectors 

For MODE^O, each filtered codevector, fj(n), is decorrelated to the long term predictor vector, b"L(n). If MODE=0, the 
t>L( n ) vector and the single VSELP codebook excitation are replaced by two VSELP codebook excitations, as the 
excitation sources. In this case, decorrelation is only performed for the second VSELP codebook. In this case, the 
orthogonalization is done with respect to the filtered codevector chosen from the first VSELP codebook. 

Defining: 

N s -\ 

T= YJb" L {n)) (110) 

n = 

and 

N s -\ 
¥m= l>" L (n)? m (rc) ( m ) 

n = 

for 1 < m < M; then q' m (n), the decorrelated filtered basis vectors, can be computed by: 

A i 



q' m (n)=q m (n)- 



in 



b" L (n) (112) 



v x J 



r 

for 1 < m < M and < n < N s -1. 

The decorrelated filtered code vectors can now be expressed as: 

M 
fi(n)=Y J im <l'm(n) (113) 

m=\ 
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for < i < 2 M -1 and < n < N s -1. 

4.1 .1 0.2 Fast search technique 

The codebook search procedure should find the codeword i which minimizes: 
Ns-1 

E' i = Y J {pM-rf' i M) 2 (114) 

n=0 

Defining : 

;v,-i 
Q = J]//(n)/>(n) (115) 

n=0 

and 

N.-l 



G t = 2l U'M 2 < 116 ) 



«=0 
then the best code vector is the one which maximizes: 

(cA 2 



G, 



(117) 



and the corresponding optimal gain is given by: 



r t = -± (us) 



I 



The search process needs to evaluate equation (117) for each codevector. The codevector which maximizes equation 
(1 17) is then chosen. Using properties of the VSELP codebook construction, the computations required for computing 
Cj and Gj can be greatly simplified. 

defining: 

R m = 2 T,Q'm(n)p(n) (119) 

n=0 

for 1 < m < M and 

N s -1 
D mj =4Y,q' m (n)q'j(n) (120) 

n=0 

for 1 <m<j<M 

Cj can be expressed as: 

1 M 

Ci=-Y< e im R m (121) 

m=\ 

and Gj can be expressed as: 

M j'-l i M 



1 iVl j — i 1 IV1 

G ; =-E E^^ wj +tI^ (122) 



7=2 m=l >=1 
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Assuming that codeword u differs from codeword i in only one bit position, say position v such that 
UV = -0i v and um = 9j m for m £ v then: 

C u = Q + 0uv R v (123) 

and 

v-1 M 

G u = G i + Y J u j0uvD jv + Ys e uAvDvj (124) 

J-l J=v+1 

The codebook search is structured such that each successive codeword evaluated differs from the previous codeword in 
only one bit position, then equation (123) and equation (124) can be used to update Q and G[ in a very efficient 
manner. Sequencing of the codewords in this manner is accomplished using a binary Gray code. This updating 
operation for both equation (123) and equation (124) can be accomplished with a total of only M multiply-accumulates 
per codevector. 

With this technique for computing Cj and Gj for the codevectors in a VSELP codebook, it is necessary to find the i 
which maximizes equation (117). Note that complementary codewords (see subclause 3.1.10) will have equivalent 
values for equation (117). Therefore only half of the codevectors need to be evaluated. Once the codevector which 
maximizes equation (1 17) is found, the sign of Cj for that codevector will determine whether that codevector or its 

complement will yield a positive gain, y 

A running maximum for equation (1 17) is kept during the code search then for each codevector evaluated, evaluate 
equation (117) and compare to the running maximum. 

-> y ' (125) 



"i '-'best 

Evaluating equation (125) directly from Cj and Gj requires one multiply, one divide and one compare operation. By 
cross multiplying equation (117) can be expressed as: 

(C i ) 2 G hest >(C hest ) 2 G i (126) 

Using equation (126) requires only three multiplies and a compare per evaluation (and no divides) where (C^est) an d 
Gbest are updated throughout the search to reflect the running best codeword. 



4.1 .1 1 Multimode gain vector quantization 



A separate GSP0 vector quantizer is derived for each of the four voicing modes. Once the frame voicing mode is 
selected, the vector quantizer, corresponding to that mode, is searched to select the excitation gains at each subframe of 
the frame. 

Although the interpretation of what the excitation sources are differs between MODE=0 and the remaining MODE 
values, the procedure for searching the gain vector quantizer is identical. In each case, the P0 term specifies the relative 
contribution of the first of the two excitation vectors to the total excitation energy at the subframe, where the first 
excitation vector is the long term prediction vector for MODE=l, 2 or 3, while the vector selected from the first of the 
two VSELP codebooks is used in the MODE=0 case. 

4.1.11.1 Coding GS and P0 

Define ex(n) to be the excitation function at a given subframe. For MODE=l, 2 or 3, ex(n) is a linear combination of the 
pitch prediction vector scaled by P, the long term predictor coefficient, and of the codevector scaled by y its gain. In 
equation form 

ex(n) = / 5c (n) + 7c 1 (n) 0<n<N s -l (127) 

where for MODE^O 



ETSI 



3GPP TS 46.020 version 7.0.0 Release 7 



34 



ETSI TS 146 020 V7.0.0 (2007-06) 



CQ(n) is the unweighted long term prediction vector, t>L(n) 
ci(n) is the unweighted codevector selected, uj(n) 

and for MODE=0 

cg(n) is the unweighted codevector selected from the first VSELP codebook, uj j(n) 

ci(n) is the unweighted codevector selected from the second VSELP codebook, ujj 2(n) 

The variable c'j(n) is a weighted version of c;(n). The power in each excitation vector is given by 

N,-l 
R x (k)=^c 2 k (n) 0<k<l (128) 

n=0 



Let R be the total power in the coder subframe excitation 
R = j3 2 R x (Q) + r 2 R x (l) (129) 

P0, the power contribution of the pitch prediction vector as a fraction of the total excitation power at a subframe, 
jB 2 R x (0) 



P0-- 



R 



where < P0 < 1 



(130) 



Define R' q (0) to be the quantized value of R(0) to be used for the current subframe and Rq(0) to be the quantized value 
ofR(0).Then: 



R'q(0) - Rq(0)previous frame 


for subframe 1 


(Ola) 


R'q(0) = Rq(0)current frame 


for subframes 2, 3, 4 


(131b) 


Let RS be 






N„ 







RS = N s R q (0)Y[(l-r 2 ) 



(132) 



i=\ 



The term GS is the energy tweak parameter defined as 

R 



R = GSRS ->GS = 



RS 



(133) 



P0 represents the fraction of the total subframe excitation energy which is due to the first codebook vector, and GS, the 
energy tweak factor which bridges the gap between R, the actual energy in the coder excitation, and RS, its estimated 
value. 

The gain bias factor %, formulated to force a better energy match between p(n) and the weighted synthetic excitation, is 
given below where. 



j=min 



2,0, 



max 



1,0, 



R 



pp 



' PoptRjw) + rlptRcAn) + 2/3 op t 7o P tRcc{o,i) 



(134) 



The weighted error equation is 

E = X 1 R pp ~ a^JGS PQ -b^GS{\- PQ) + cGS^P^l- PQ) + d GS PQ + eGS{\- PQ) (135) 

where 
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'' = 2 *M°),r/^ d36) 



b = 2 Z R pc (l)l-^- (137) 



'^(1) 
2R CC {0,1)RS 

V^(o)^(i) 

RSR cc {0,0) 



(138) 



(139) 



e = °£±J-L (140) 



**(o) 

A^-l 
* pc (*) = JpWc'iW k=0,l (141) 



n=0 



AT,-1 



^c(^i)=S^(n)c'jW k=0,l,j=k,l (142) 

R cc (k,j) = R cc (j,k) (143) 

v = Z^ 2(n) (144) 

Four separate vector quantizers for jointly coding P0 and GS are defined, one for each of the four voicing modes. The 
first step in quantizing of P0 and GS consists of calculating the parameters required by the error equation: 

R cc (kj) k = 0, 1, j=k, 1 

R x (k) k = 0, 1 

RS 

Rp C (k) k = 0, 1 

a, b, c, d, e 

Next equation (135) is evaluated for each of the 32 vectors in the {P0,GS} codebook, corresponding to the selected 
voicing mode, and the vector which minimizes the weighted error is chosen. Note that in conducting the code search 

J^Rpp may be ignored in equation (135), since it is a constant, pq, the quantized long term predictor coefficient, and 

Yq, the quantized gain, are reconstructed from 



RSGS vq P0 vq 

4=V *,(o) (145) 



\RSGS vq \l-P0 vq , 

-^ (146) 

where P0 V q and GS vq are the elements of the vector chosen from the {P0,GS } codebook. 
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A special case occurs when the long term predictor is disabled for a certain subframe, but voicing MODE^O. This will 
occur when the state of the long term predictor is populated entirely by zeroes. 



For that case, the following error equation is used: 

E^^Rpp-b^GS+eGS 

For this case the quantized codevector gains are: 

Pq=0 



(147) 



(148) 



\RSGS vq [l-P0 vq 



(149) 



4.2 GSM half rate speech decoder 
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Figure 5: The GSM half rate speech decoder for MODE = 1, 2 or 3 

A block diagram of the GSM half rate speech decoder for MODE=l, 2 or 3 is given in figure 5. The speech decoder 
creates the combined excitation signal, ex(n), from the long term filter state and the VSELP codevector. For MODE=0, 
the long term filter state is replaced by another VSELP codebook and the pitch prefilter is not used. The combined 
excitation is then processed by an adaptive pitch prefilter and gain. The prefiltered excitation is applied to the LPC 
synthesis filter. After reconstructing the speech signal with the synthesis filter, an adaptive spectral postfilter is applied 
followed by an automatic gain control which is the final processing step in the speech decoder. 
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4.2.1 Excitation generation 

The combined excitation, ex(n), shall be computed as shown in equation (127) 

The combined excitation, ex(n), is filtered by the synthesis filter to generate the speech signal. The synthesis filter is a 
tenth order all pole filter. The filter coefficients for the subframe are the aj's defined in subclause 4.1.6. The filter 
coefficients will change from subframe to subframe. The filter state shall be preserved from subframe to subframe. A 
direct form filter shall be used for the synthesis filter. 

4.2.2 Adaptive pitch prefilter 

Given ex(n) as the input, exp(n), the pitch prefiltered output, is defined by 

ex p (n) = ex(n) + ^ex p (n-L) ;for0<n<N s -l (150) 

where 

, Jo,3Mj«[/?,VPOJ ; MODE * 
I ;MODE = 

Since L can be fractional in value, an interpolating filter is used. This is the same interpolating filter which is used for 
the open loop lag search. A gain scale factor is computed and is used to scale the pitch prefiltered excitation, prior to 
applying it to the LPC synthesis filter. Pscale* me gain scale factor, is 



"scale 



N-l 



J\x 2 (w) 



^ (152 ) 



n=0 
Thus exp S (n), the gain corrected pitch prefiltered excitation which drives the LPC synthesis filter, is given by 

ex pA n ) = P scale ex p{ n ) ; for < n < N s -1 (153) 

4.2.3 Synthesis Filter 

A direct form synthesis filter is used: 
10 

s(n) = ex ps (n) + ^ j a i s(n-i) ,0<n<N s -l (154) 

(=1 



4.2.4 Adaptive spectral postfilter 

The perceptual quality of the synthetic speech is enhanced by using an adaptive postfilter as the final processing step. 
The general form of the postfilter is given by: 

N P . 
s{n) = s{n)-Y.rfcc i s(n-i) ,0<n<N s -l (155) 

i=\ 



l(n) = s(n) + ^ i {0,75)'a i s(n-i) ,0<n<N s -l (156) 



/=i 
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s(ri) = s(ri)-0,2s(n-l) 



0<n<N s -l 



(157) 



The adaptive spectral postfilter numerator polynomial equation (155) is replaced by a spectrally smoothed version of 
the adaptive spectral postfilter denominator polynomial equation (156). To derive the coefficients of the numerator 
polynomial, the denominator polynomial coefficients are converted to the autocorrelation coefficients R(i). The SST 
bandwidth expansion function is then applied to the autocorrelation sequence, 



R sst (i) = R(i)w sst (i) 



, < i < N r 



(158) 



and the numerator polynomial coefficients are calculated from the modified autocorrelation sequence via the AFLAT 
recursion. 

From R ss t(i) the reflection coefficients which define the combined spectrally noise weighted synthesis filter are 
computed using the AFLAT recursion once per frame. 



STEP1 



Define the initial conditions for the AFLAT recursion: 
P (i) = R sst (i) ,0<i<N p 



V„d) = R,Ji + l 



,l-Np<i<N p -l 



(159) 
(160) 



STEP 2 Initialize j, the index of the lattice stage, to point to the first lattice stage: 



STEP 3 



j=l 

Compute r;, the j-th reflection coefficient, using: 

V)-l(0) 



(161) 



r J = 



r J-l 



(0) 



(162) 



STEP 4 



.0<i<Np-j-l (163) 



STEP 5 

STEP 6 
STEP 7 



Given r;, update the values of Vj and P; arrays using: 

Pj (0 = (l + r})Pj_! (i) + rj \Vj_i (i) + V M (-») 

Vj (0 = Vj-! ( i + 1) + rjVj_! (- i - 1) + 2 rj Pj_ { (| i + 1|) , 1+j -Np<i<Np-j - 1 ( 1 64) 

Increment j : 

j=j+l 

If j < Np go to step 3, otherwise all Np reflection coefficients have been obtained. 

The reflection coefficients, rj, are then converted to OC v the direct-form LP filter coefficients for use in 
the adaptive spectral postfilter numerator polynomial. 



The resultant adaptive spectral postfilter is derived from equations 155, 156 and 157: 



s(n) = s(n)— Y.OiiS{n-i) 
i=\ 

N p 
s(n) = s(n)+ I (0,75)* a t s(n-i) 
i=\ 

s(ri) = s(n)-0,2s(n-l) 



,0 < n < N s -1 



(165) 



,0 < n < N„-l 



(166) 



,0<n<N s -l 



(167) 



In order to reduce the computations needed to compute the spectrally smoothed numerator coefficients, the spectral 
smoothing operation is performed once per frame on the denominator coefficients corresponding to the uninterpolated 
coefficients. This will yield the coefficients for the numerator of the spectral postfilter for sub frame four. The numerator 
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coefficients for subframes one, two, and three are interpolated using the same interpolation scheme that is used for the 
LPC synthesis coefficients (see subclause 4.1.6). 

As in the case of the pitch prefilter, a means of automatic gain control is needed to ensure unity gain through the 
spectral postfilter. A scale factor, S sca le is given by: 

S'scale(n) = (0,9875 S' sca i e (n-1) ) + (0,0125 S sca i e ) (168) 

Scale factor, S sca j e is the square root of the ratio of the input signal energy to the output signal energy over the 
subframe 

The output of the spectral postfilter is then multiplied by S' sca i e as the last step in reconstructing the speech signal in the 
speech decoder. 

4.2.5 Updating decoder states 

The long term predictor state, r(n), is updated by: 

r(n) = r(n+40) for -146 < n <-41 (169) 

r(n) = ex(n+40) for -40 < n < - 1 (170) 



5 Homing sequences 

5.1 Functional description 

The half rate speech codec as well as the DTX system and comfort noise generator are described in a bit exact 
arithmetic to allow for easy type approval as well as general testing purposes of the half rate speech codec. 

The response of the codec to a predefined input sequence can only be foreseen if the internal state variables of the codec 
are in a predefined state at the beginning of the experiment. Therefore, the codec has to be put in a so called home state 
before a bit exact test can be performed. This is usually done by a reset. 

To allow a reset of the codec in remote locations, special homing frames have been defined for the encoder and the 
decoder, thus enabling a codec homing by inband signalling. 

The codec homing procedure is defined in such a way, that on either direction (encoder or decoder), the homing 
functions are called after processing the homing frame that is input. The output corresponding to the first homing frame 
is therefore dependent on the codec state when receiving that frame and hence usually not known. The response to any 
further homing frame in one direction is by definition a homing frame of the other direction. This procedure allows 
homing of both, the encoder and decoder from either side, if a loop back configuration is implemented, taking proper 
framing into account. 

5.2 Definitions 

encoder homing frame: The encoder homing frame consists of 160 identical samples, each 13 bit long, with the least 
significant bit set to "one" and all other bits set to "zero". When written to 16 bit long words with left justifications, the 
samples have a value of 0008 hex. Test sequence SEQ05.INP described in GSM 06.07 [3] defines the encoder homing 
frame. The speech decoder has to produce this frame as a response to the second and any further decoder homing frame 
if at least two decoder homing frames were input to the decoder consecutively. 

decoder homing frame: The decoder homing frame has a fixed set of speech parameters as defined in test sequence 
SEQ05.INP described in GSM 06.07 [3], It is the natural response of the speech encoder to the second and any further 
encoder homing frame if at least two encoder homing frames were input to the encoder consecutively. 
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5.3 Encoder homing 



Whenever the half rate speech encoder receives at its input an encoder homing frame exactly aligned with its internal 
speech frame segmentation, the following events take place: 

Step 1: The speech encoder performs its normal operation including VAD and DTX and produces a 

speech parameter frame at its output which is in general unknown. But if the speech encoder was 
in its home state at the beginning of that frame, then the resulting speech parameter frame is 
identical to the decoder homing frame (this is the way how the decoder homing frame was 
constructed). 

Step 2: After successful termination of that operation, the speech encoder provokes the homing functions 

for all submodules including VAD and DTX and sets all state variables into their home state. On 
the reception of the next input frame, the speech encoder will start from its home state. 

NOTE: Applying a sequence of N encoder homing frames will cause at least N-l decoder homing frames at the 
output of the speech encoder. 

5.4 Decoder homing 

Whenever the speech decoder receives at its input a decoder homing frame, then the following events take place: 

Step 1: The speech decoder performs its normal operation including comfort noise generation and 

produces a speech frame at its output which is in general unknown. But if the speech decoder was 
in its home state at the beginning of that frame, then the resulting speech frame is replaced by the 
encoder homing frame. This would not naturally be the case but is forced by this definition here. 

Step 2: After successful termination of that operation, the speech decoder provokes the homing functions 

for all submodules including the comfort noise generator and sets all state variables into their 
home state. On the reception of the next input frame, the speech decoder will start from its home 

state. 

NOTE 1: Applying a sequence of N decoder homing frames will cause at least N-l encoder homing frames at the 
output of the speech decoder. 

NOTE 2: By definition the first 58 bits of the decoder homing frame must differ in at least one bit position from the 
first 58 bits of any of the decoder test sequences. Therefore, if the decoder is in its home state, it is 
sufficient to check only these first 58 bits to detect a subsequent decoder homing frame. This definition 
ids made to support a delay optimised implementation in the TRAU uplink direction. 



5.5 Encoder home state 



In GSM 06.06 [2], a listing of all the encoder state variables with their predefined values when in the home state is 
given. 



5.6 Decoder home state 



In GSM 06.06 [2], a listing of all the decoder state variables with their predefined values when in the home state is 
given. 
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Annex A (normative): 
Codec parameter description 

A. 1 Codec parameter description 

The following is a list of all the parameters which are coded for each 20 ms speech frame. The basic data rate of the 
speech coder is 5,6 kbps. Therefore each 20 ms speech frame consists of 1 12 bits. These bits are given in table A. 1 . 

Table A.1 : Codec parameter description 



Parameter | No. of bits \ Description 


Frame bits: 


MODE 


2 


voicing mode 


R0 


5 


frame energy 


LPC1 


11 


reflection coefficient vector r-j -r3 


LPC2 


9 


reflection coefficient vector r 4 -rg 


LPC3 


8 


reflection coefficient vector r 7 -r-| q 


INT LPC 


1 


the soft interpolation bit for the frame 


Subframe bits (MODE = 1 ,2 or 3): 


LAG 1 


8 


lag for first subframe 


LAG 2 


4 


lag delta code for second subframe 


LAG 3 


4 


lag delta code for third subframe 


LAG 4 


4 


lag delta code for fourth subframe 


CODE 1 


9 


codebook, I, for first subframe 


CODE 2 


9 


codebook, I, for second subframe 


CODE 3 


9 


codebook, I, for third subframe 


CODE 4 


9 


codebook, I, for fourth subframe 


GSP0 1 


5 


{P0,GS} code for first subframe 


GSP0 2 


5 


{P0,GS} code for second subframe 


GSP0 3 


5 


{P0,GS} code for third subframe 


GSP0 4 


5 


{P0,GS} code for fourth subframe 


Subframe bits (MODE=0): 


CODE1 1 


7 


codebook code, I, for first subframe 


CODE2 1 


7 


codebook code, H, for first subframe 


CODE1 2 


7 


codebook code, I, for second subframe 


CODE2 2 


7 


codebook code, H, for second subframe 


CODE1 3 


7 


codebook code, I, for third subframe 


CODE2 3 


7 


codebook code, H for third subframe 


CODE1 4 


7 


codebook code, I, for fourth subframe 


CODE2 4 


7 


codebook code, H, for fourth subframe 


GSP0 1 


5 


{P0,GS} code for first subframe 


GSP0 2 


5 


{P0,GS} code for second subframe 


GSP0 3 


5 


{P0,GS} code for third subframe 


GSP0 4 


5 


{P0,GS} code for fourth subframe 



A.1.1 MODE 

The speech coder is defined by 4 voicing modes. MODE is a two bit code which specifies which of the four voicing 
modes is used at the current frame. The MODE indicates which definition of the frame bits to apply to the current 
frame. 

A.1.2 R0 

R0 is a code which represents the average signal power of the input speech for the frame. The average signal power is 
computed using an analysis window which is centered over the last 100 samples of the frame. 
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A.1.3 LPC1-LPC3 

The 10 reflection coefficients are vector quantized in three vector segments. The first vector segment codes reflection 
coefficients i\ - ij, the second vector segment codes coefficients rz|.-r6, the third vector segment codes coefficients r7 - 

riO- 

A. 1.4 LAG_1 -LAG_4 

LAG_1, the lag for the first subframe, can take on the value in the range of 21 to 142. Eight bits are used to encode the 
lag which may be fractional in value. Each of the remaining lag values ( LAG_2 through LAG_4) is delta coded relative 
to the preceding subframe's coded value of the lag, with a deviation of -8 to +7 allowable lag value levels specified by a 
four bit code. 

A.1 .5 CODExJ - CODEx_4 

If MODE ± 0, the code value for the VSELP codebook is the codeword I as derived by the codebook search procedure. 
If MODE=0, two VSELP codebooks are sequentially searched, with codeword I, specifying the codevector from the 
first VSELP codebook, assigned onto CODEl_x, and codeword H, specifying the codeword selected from the second 
VSELP codebook, assigned onto CODE2_x, where x is the subframe number. 

A.1.6 GSPOJ -GSP0_4 

The { P0,GS } codebook contains the values needed to determine the gain factors for the excitation vectors of a given 
subframe. The index of the corresponding codebook entry is assigned to GSP0_x. 

The speech coder is a multimode speech coder, defined by four voicing modes: 

MODE = unvoiced 

MODE = 1 slightly voiced 

MODE = 2 moderately voiced 

MODE = 3 strongly voiced 

If MODE=0, the adaptive codebook (long-term predictor) and the VSELP codebook are replaced by two other VSELP 
codebooks. 



A.2 Basic coder parameters 

The following are the basic parameters for the 5 600 bps GSM half rate speech codec system. 





sampling rate 


8 kHz 


N F 


frame length 


160 samples (20 ms) 


Ns 


subframe length 


40 samples (5 ms) 


Np 


short term predictor order 


10 
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Annex B (normative): 

Order of occurrence of the codec parameters over Abis 

The order of occurrence of the codec parameters over the Abis is defined for unvoiced speech (MODE = 0) and voiced 
speech (MODE = 1, 2 or 3) in tables B.l and B.2 respectively. 

Table B.1 : Occurrence of the codec parameters over Abis for unvoiced speech (MODE = 0) 



Parameter 


No. of bits 


Bit No. (MSB - LSB) 


R0 


5 


b1 -b5 


LPC1 


11 


b6-b16 


LPC2 


9 


b17-b25 


LPC3 


8 


b26 - b33 


INT LPC 


1 


b34 


MODE 


2 


b35 - b36 




CODE1 1 


7 


b37 - b43 


CODE2 1 


7 


b44 - b50 


GSP0 1 


5 


b51 - b55 




CODE1 2 


7 


b56 - b62 


CODE2 2 


7 


b63 - b69 


GSP0 2 


5 


b70 - b74 




CODE1 3 


7 


b75 - b81 


CODE2 3 


7 


b82 - b88 


GSP0 3 


5 


b89 - b93 




CODE1 4 


7 


b94-b100 


CODE2 4 


7 


b101 -b107 


GSP0 4 


5 


b108-b112 



Table B.2: Occurrence of the codec parameters over Abis for voiced speech (MODE = 1, 2 or 3) 



Parameter 


No. of bits 


Bit No. (MSB - LSB) 


R0 


5 


b1 -b5 


LPC1 


11 


b6-b16 


LPC2 


9 


b17-b25 


LPC3 


8 


b26 - b33 


INT LPC 


1 


b34 


MODE 


2 


b35 - b36 




LAG 1 


8 


b37 - b44 


CODE1 


9 


b45 - b53 


GSP0 1 


5 


b54 - b58 




LAG 2 


4 


b59 - b62 


CODE2 


9 


b63 - b71 


GSP0 2 


5 


b72 - b76 




LAG 3 


4 


b77 - b80 


CODE3 


9 


b81 - b89 


GSP0 3 


5 


b90 - b94 




LAG 4 


4 


b95 - b98 


CODE4 


9 


b99-b107 


GSP0 4 


5 


b108-b112 
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